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ABSTRACT 

Mixed modes seen in evolved stars carry information on their deeper layers that can place stringent 
constraints on their physics and on their global properties (mass, age, etc.). In this study, we present 
a method to identify and measure all oscillatory mode characteristics (frequency, height, width). 
Analyzing four subgiants stars, we present the first measure of the effect of the degree of mixture on 
the £ = 1 mixed modes characteristics. We also show that some stars have measurable £ = 2 mixed 
modes and discuss the interest of their measure to constrain the deeper layers of stars. 
Subject headings: stars: oscillations, stars: interiors, methods: data analysis 



1. INTRODUCTION 

Identifying oscillatory modes and measuring their 
properties is fundamental to asteroseismology. 
For solar-type stars on the main-sequence, the observed 
oscillations are pressure waves (or p modes), while grav- 
ity modes (or g modes) are confined to the radiative 
core and have a surface amplitude too small to be mea- 
sured directly (jBelkacem et al.|[2009l : lAppourchaux et ahl 
l201Clf ). Therefore, the frequency spectrum follows the 
asymptotic rela tion for the p modes which is a rather 
simple pattern ()Vandakurovl(l967l : iTassoullHOSOl: iGoughl 
[T98@ . 

As a star evolves into a subgiant, the deeper layers of 
the core stop burning hydrogen which modifies signifi- 
cantly the internal structure and therefore the stellar os- 
cillations. Modes that exhibit both p and g modes p rop- 
erties, so-called mixed modes (jAizenman et aLlflQ??! ). be- 
come visible in power spectra of such stars. 

These mixed modes can be hard identify and analyze 
because they do not follow the asymptotic relations of 
pressure-modes. However, a recent dedicated asymptotic 
development has been proposed bv lMosser et al.l (|2012cD 
easing the mode identification. 

Mixed modes are sensitive to the deeper layers of stars 
and allow us to probe more precisely their evolution- 
ary state. For instance, rapid changes in the core of 
evolved stars modify the gravity-mode frequencies in 
very short timescales. Consequently, the fr equencies of 
the mixed modes chan ge quickly with time ()Osakilll975t 
lAizenman et al.lll977[ ). potentially providing stell ar ages 
with a precision down to a few Myr (Metcal fe et al.l 
[20T0I) . 

Until now, the measures of the individual modes 
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characteristics were based on assumptions valid i n the 
Sun and similar stars (e .g. lAppourchaux et ahl 120081 
lAppourchaux et al.ll2012| ). Although these assumptions 
are well suited for low-mass main-sequence stars, they 
do not hold for more evolved stars. For example, an 
in-dep th theoretical study of red giants by iDupret et al.l 
(|2009f ) showed that widths and amplitudes of the non- 
radial modes are smaller than for radial modes and may 
vary greatly from one mode t o another as the modes un- 
dergo avoided crossing s {e.g. IDeheuvels fc MichelllMll : 
IBenomar et al.|[2"012b| ). A precise measure of the prop- 
erties of the non-radial modes would therefore provide 
unprecedented constraints on the damping and on the 
inertia of modes confined to the deeper layers of stars. 

The goal of this paper is to describe a global method 
identifying the non-radial modes and measuring all indi- 
vidual mode properties in subgiants observed by Kepler^ 
using an appropriate set of assumptions. A discussion 
with theoretical expectations is also presented. 

We first describe the potential of mixed modes as 
probes of the stellar properties (Section We then 
develop the methods and formalism for measuring the 
properties of the mixed modes (Section [3|), define the 
limitations of the approach (Section and finish by pre- 
senting (Section [5|) and discussing (Section (5]) the results 
from observations. 

2. MIXED MODES IN EVOLVED STARS 

Mixed modes have an accoustic behavior in the 
envelope of the star and behave as gravity waves 
in the deep interior; the relative importance of 
the acoustic and gravity-wave behavior is deter- 
mined by the coupling in the mode between these 
two regions, which depends on the frequency and 
degree of the mode. They occur in evolved stars (sub- 
giants and red giants), in which the large density gradi- 
ent outside the core effectively divides the star into two 
coupled cavities. Mixed modes have p-mode character 
in the envelope and g-mode character in the core. This 
leads to mode bumping, in which mode frequencies are 
shifted from their regular spacing and no longer follow 
the asymptotic relations for p modes that is seen in main- 
sequence stars. 

We can model mixed modes as arising from fictitious 
pure p and g modes that would exist if their cavities 
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were not coupled. In o rder to avoid ambiguity we follow 
lAizenman et alj ()1977f ) by referring to t hese pure mod es 
as TT and 7 modes respectively (see also lBedding||20lH ). 

Mode bumping in subgiant stars was first ob- 
serve d and modeled in rj Boo feieldsen et alJ 
19951 iChristensen-Dalseaard et al.l 



Guenther fc Demargud 119961 : iKieldsen et alJ 
20051) and /3 Hyi (iBeddiilg et alJ 



Carrier et al 



Brandao et al 



1995 



2003 



2007 



|2011|) . More recently, asteroseismic space 



missions have produced many more exa mples, inclu ding 
the CoRoT target HD 49835 (fPcheuvcls et al.li2010D and 



20101: iMathur et alJ 


2OIOI: iCampante et al.l 
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: lAppourchaux et al.l 



Meanwhile, thousands of giants have been observed with 
CoRoT and Kepler showing a 'forest' of non-radial modes 
(IBedding et al.llMol: iHuber et al.|[2010l : IKallinger etHI 
I20TI " identified as = 1 mixed modes. This has led to 
two major discoveries: a seismic diagnostic that permits 
to distin guish RGB stars from red clump stars 



us 



(iBeck et al.l I2OIII: IBedding et all 120111: iMosser et ahl 

1201 1[) and e vidence for radial differential rotation 

within stars (IBeck et all l2012t IPeiieuvels et al, 20l1 
Mosser et all l2012dl : iMarques et al.l 120131: iGoupil et al l 



Although all non-radial modes can have mixed prop- 
erties, most work to date on subgiants and giants has 
focused on dipole modes = 1) because the distance 
between the pressure and gravity wave cavities is less 
than for higher-degree modes. Hence the trapping 
is less efficient for dipole modes than at higher 
degrees. Observationally it means that dipole 
modes l ifetimes are smaller and their amplitudes 
higher (jPupret et al.l l2009l ): therefore they are 
more easy to resolve and to observe. 

In the early subgiant phase, only few avoided cross- 
ings are observed because the 7 modes are well spaced 
in frequency (low 7-mode density). As a subgiant 
evolves, the 7-modc density increase, thus the num- 
ber of 7 modes within the observable range of fre- 
quency increases (approximately between 0.5t'inax 
and 1.5i^inax5 where i^max is the frequency at max- 
imum power for the modes). This leads to numer- 
ous (measurable) avoided crossings. With more avoided 
crossings, the spectrum becomes very complicated, which 
makes mode identification difficult, because mixed modes 
spread all over the echelle diagrams. When the number 
of TT and 7 modes become similar, the spectrum is par- 
ticularly difficult to interpret. 

Recently, theoretical developmen ts (iGoupill 120131 i n 
prep.) following the formalism of lUnno et al.l ( 19891 ) . 
have led to an asymptotic relation for mixed modes, 
which has been applied by iMosser et al.l (j2012cl ) to gi- 
ant stars. With this asymptotic expansion, modes can 
be more easily identified in giants and subgiants. The 
asymptotic relation for mixed modes uses constant val- 
ues of the large separation (Aiy) and of the period spac- 
ing (Allf), to characterize the mixed modes and is well 
suited to measuring global frequency characteristics of 
the mode bumping. However it does not provide other 
mode properties (e.g. height or width of the modes). 

The approach presented in this paper aims to fit all 
individual properties of the modes. It describes the fre- 
quencies of £ = 1 modes as eigensolutions of a series 



of coupled harmonic oscillators, w hich has been shown 
to fit observed avoided cro ssings (jPeheuvels fc Michell 
I2OIII : iBenomar et al.l l2012al) . The approach uses a full 
Bayesian framework to include our prior knowledge of 
the pressure and gravity modes. 

3. FITTING THE POWER SPECTRUM 

In this section, we present the model that describes the 
frequencies, referred as the coupled-oscillators model and 
the method and assumptions for the fit of the spectra. 

3.1. The Coupled- Oscillators Model 

In the coupled oscillator model we expand 
the oscillations of a stellar model on fictitious 
pure TT and 7 modes, with frequencies j = 

{1, . . . , iVTr}? and o;^'*''', k = {I, . . . , N^}, respectively. 
We seek angular frequencies uj = {ui, . . . ,ujn}, 
with N = + N^, which are to be fitted to the 
observed frequencies. This leads to the following 
system of equation, 

AY = uj^Y, (1) 

where Y = {yi,y2, ■■■jUn} is a vector of mode ampli- 
tudes and A is a matrix that contains the frequencies of 
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(2) 

Note that the coupled oscillator does not repre- 
sent a physical model of stellar oscillations but 
it does contain the salient aspects of the prop- 
erties of the mixed modes, allowing it to pro- 
vide a description of the observed frequencies in 
terms of the relevant pr o perti es of the star. Fol- 
lowing [Peheuvels fc Micheil ()2010f ) and IBenomar et al.l 
(l20T2bl) . we assumed for each 7 mode (that is, at 
each avoided crossing), that the coupling strength be- 
tween this 7 mode and the iV^r tt modes has a sin- 
gle value, cki, with i = {1,...,A^^}. Note that at 
early stages of the sub-giant phase, TV^ ^ A^-y: the 
spectrum contains few avoi ded crossings and has a 
p-mode dominated pattern (iDeheuvels fc Michell 120111 : 
ICampante et al.ll201l"l: lAppourchaux et al.ll2012D . At the 
other extreme, red giants have ^ and the power 
spectrum has a g-mode do minated pattern (iDupret et al. 
20091: IBedding 2011: StcUo 20 lit iChristensen-Dalsgaaro 
2012l : [Mosser et al.,,2012G ). Between these two extremes. 



one can consider inte rmediate cases where A^^r ~ (e.g., 
Idi Mauro eraI1l201lD . whose power spectra may be hard- 
est to interpret. 

3.2. Markov Chain Monte Carlo fitting 

In recent years, Markov Chain Monte Carlo (MCMC) 
approache s have been widely used in asteroseismol- 
ogy (e.g. iBrewer et all 120071: iGruberbauer et al.l 120091: 
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Benomar etldl |2009allbl: iHandberg fc Campanti 120111: 



Campante et alj 120111 : lAppourchaux et al.ll2012[ ). Such 
approaches rely on samphng the posterior Probabihty 
Density Function (PDF) and provide far more robust re- 
sults than maximization methods such as the Maximum 
Likelihood Estimator (MLE) or the Maximum A Poste- 
riori (MAP) approach. 

MCMC fitting is well suited to the present analysis be- 
cause it provides an easy way to compute the PDF of any 
function of the fitted variables. For example, the PDF 
of the mean large separation or the mean period spacing 
(cf. Fig. [T]) can be easily derived from the samples ac- 
quired by MCMC, for each individual p-mode or g-mode 
frequencies. When dealing with very complex functions 
or highly non-linear models, such as the coupled oscilla- 
tors model, it provides an easy way to propagate the un- 
certainties. Note that the matrix A contains often hun- 
dreds of correlated parameters, so the error propagation 
law would require us to determine thousands of deriva- 
tive terms, either analytically or numerically. MCMCs 
allow one to propagate the error in a more straightfor- 
ward and simpler way. Moreover it does not rely on the 
PDFs being Gaussian. 

3.3. Defining the posterior probability density function 

In a Bayesian approach, the first step is to write the 
likelihood function and define the posterior PDF. We 
deal with a power spectrum, thus the noise statistic of 
each data point yi follows a with two degrees of free- 
dom. Knowing the statistics of data points is sufficient 
to derive the likelihood function, give n the fitted model 
M(iy, 9) that depends on variables 6 (jPuvall fc HarvevI 
[1981 : 



N 



p{y\e,M,I) = l[ 



1 



6) 



where TV is the number of frequencies Ui and / encom- 
passes all other contextual information. Equation ([3]) is 
only applicable if the sampled frequencies are indepen- 
dent, that is, if the power spectrum is sampled at the 
formal frequency resolution of the data set. 

From Bayes theorem, we can define the posterior PDF 
that we seek to sample. 



p{e\y, M, /) = p{e\M, I)p{y\e, M, I)/C 



(4) 



where p{0\M,I) is our explicit prior knowledge on the 
parameter set 6 of the model A/, and C is a normalization 
constant. To apply priors on a particular parameter (or 
a subset of parameters), we use the product rule for the 
case of independent variables: the global prior is simply 
the product of the individual priors. For example, we 
can separate the prior associated with the noise On from 
the prior on the mode 0s, 

p{es,ON\Mj)^pies\Mj)pieN\M,i). (5) 

3.4. Defining the priors 
3.4.1. Frequencies 

To a first approximation, each resolved mode is well 
represented by a Lorentzian in the power spectrum. This 
Lorentzian is characterized by a height H, a width F 



and a central frequency u. The frequency distribu- 
tion of TT modes is approximately described by 
the sec ond-orde r asymptotic relation for pressure 
modes ()TassouH [l98'0f), 



£ 



- Si 



(6) 



where Si^oi represents the second-order term and 
~ 1/4 fo r low-m ass stars, as verified by 
osser et al.l ()2012al ). For an accurate descrip- 
tion, one must take care to distinguish the ob- 
served and asymptotic values of the large separa- 
tion. For 7 modes, the asymptotic relation for periods 
P-v is (|Tassoullll98(l 



) An., 



a-y) Allf, (radiative core) 



(convectivc core). 



(7) 

Here and are the radial orders for tt modes and 7 
modes, (. is the degree, is a phase offset related to the 
position of the waves turning point in the stellar atmo- 
sphere (e.g. IWhite et al.]l2011[ ). while a-y and 5vm carry 
information about the core. AIX^ is the period spacing of 
7 modes with degree £ such as Allf = Alio / -\/«(Z+l) . 

These equations are used in order to define the priors 
on the frequencies of the tt and 7 modes. Note that, in 
practice, cannot be directly measured as we can't dis- 
tinguish a priori between convective and radiative cores. 
Therefore, we will use e^^i to denote the phase offset for 
the gravity-modes, independent of the nature of the core. 

Abrupt changes in the physical properties inside stars, 
known as glitches, cause departures from the predicted 
asymptotic frequencies for both pressure and gravity 
modes. We ex pect the depar tures to be less than a few 
percent of At^ ilProvost et al.|[T993HMazumdar fc Michell 



(3) [2Mnl : IMiglio et al.ll201fln . 



In order to take these glitches into account (and to 
reproduce accurately the observed mixe d mode frequen- 
cies) , w e used an approach similar to iBenomar et al.l 
(|2012a[) for the tt modes. We applied smoothness con- 
ditions on the average large separation, as well as on the 
second derivative in tItt for the frequencies iy^(n7r). The 
priors on Av ensure that the £ = 1 tt modes stay approx- 
imately aligned in the echelle diagram, while the prior on 
the second derivative avoids strong, non-physical kinks. 
This was achieved by applying a Gaussian prior on these 
quantities. For example, for the second derivative in fre- 
quency A'^v{n), we have, 



p{A^v{n)\M,I) 



■ exp 



a 1^1 



(8) 

where ai^iy plays the role of a relaxation constraint 
and must be chosen to ensure enough freedom, but not 
too much, in order to efficiently smooth the frequency 
profile. Here, a trial-and-error procedure showed that 
values of (TA^iy between 0.5 to 1.5 /iHz offer a good com- 
promise: stars with a greater number of mixed modes 
need a tighter prior. 

In a very similar way, smoothness conditions were ap- 
plied to the 7 modes, with the difference that they were 
applied to the periods of the modes (period spacings and 
second derivative of the periods), instead of frequencies. 
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In subgiants the number of 7 modes within the range of 
observed frequencies is small (i.e. N^^ > N^). There- 
fore to avoid overfitting, the relaxation constraint was 
set such that mode-to-mode variations of more than 1% 
in period spacing were unlikely. 

With such smoothness conditions, the 7r-mode and 7- 
mode deviation from a strictly regular pattern of fre- 
quencies or periods is locally described by a second-order 
polynomial function of the radial order, and the solution 
belongs to the family of spline functions. 

3.4.2. Other parameters 

We also need to define priors for the mode linewidths. 
Mode bumping not only affects modal frequencies, but 
also their widths. Indeed, the higher the mode ampli- 
tude in the 7-cavity , the smaller the mode linewidth (e.g. 
IDupret et al.ll2009f i. Moreover, mixed modes close to the 
i = 1 TT-mode 'home ridge' {i.e. the expected frequen- 
cies of pure p modes), have more pressure-like properties 
than more strongly bumped mixed modes, which have 
more gravity-like properties. Thus, lifetimes of mixed 
modes are expected to vary significantly, and will not 
follow the usual profile seen for pure p modes in main- 
sequence stars (a decreasing function of frequency, with 
a central plateau near t'max)- 

For dipolar modes (^ = 1), it is reasonable to use the 
adjacent £ = widths as an upper limit because the 
mixed modes should be narrower than pure p modes. 
Here, we used a smooth edge (rather than an abrupt 
one, as it would be with a uniform prior), 

r if Ti < 0, 

p{Ti) = { Cifo<ri <ro, (9) 

[ cxp-^(ri-ro)V-^ ifTi >ro, 

where C = (Fq — Fi -|- V27r(T/2)~^, ensures a proper 
normalization and a = 0.1 Fq 

For i > 2 modes, we assumed that the linewidths fol- 
low those of the £ = modes and we interpolated them. 
This is justified by the fact that pressure and gravity 
modes for £ = 2 modes are less coupled because the 
width of evanescent region between the cavi ties is larger 
than for i = 1 modes (e.g., Eq. (16.51) of lUnno et al.l 
I1989D . Moreover, in the 7 cavity the radiative damping 
increases as £^ fe.g.. iGodart et al.|[2009f ). Mixed-mode 
amplitudes decrease quickly as £ increases and modes far 
from their home ridge are difficult to observe for £ > 2. 
Consequently, for £ > 2, observed modes are likely to be 
p-dominated and exhibit mode lifetimes comparable to 
£ = modes (iDupret et al.ll2009D . 

For mode heights, which are correlated with their 
widths, we chose to simply use Jeffrey (non-informative) 
priors (it would introduce biases to use informative priors 
for both widths and heights). 

The noise background characterization can be critical 
when dealing with low signal-to-noise data. Convective 
motions at the surface of a star are in general the main 
sources of noise when we analyze the modes of solar-like 
stars in phot ometry. In th e present study, two Harvey- 
like profiles ()Harvevlfl985l ) plus white noise are used in 
order to accou nt for the different scales of the convec- 
tion cells {e.g. iMathur et al.|[20Tl[) . Priors on the noise 
parameters were based on a fit where the mode enve- 
lope was described by a Gaussian function, e.g. as in 



iBenomar etldl (|2009a[) . 

Rotation and magnetic fields are known to lift the de- 
generacy on the degree £ such that each degree shows 
a fine structure and is split into 2£ + 1 modes. In 
most cases the magnetic field is negligible, so this 
splitting is inter preted as a sig nature of stellar rota- 
tion ( e.a iGizon fc Solanki 2003 . iDeheuvels et al.l I20T2I: 
Mosser et aLn2012di: iMargues et al.l 120131 : iGoupil et al.l 
2013() . Even if the rotational splitting could be measured, 
we preferred not to implement it in this analysis. Indeed, 
several effects may make the rotational splitting hard 
to measure. First of all, it is strongly correlated with 
the inclination angle between the line of sight of the ob- 
serve r and the rotation axis of the star (jGizon fc Solankil 
I2003D . Moreover, as for the width, i t varies as func- 
tion of the frequency f or mixed modes (jBeck et al.ll2012l : 
IDeheuvels et al.ir2012f ). Finally, the height of the split 
components could be affected by the mixed properties of 
the modes. Incorporating all these possible effects would 
greatly complicate the fitted model and goes beyond the 
scope of this work. 

3.5. Defining the initial guesses 

Before fitting of the power spectrum, one needs to de- 
fine the initial guesses of the fitted model. Because mod- 
els often have hundreds of parameters, several steps are 
necessary: 

1) We use d the algorithm described bv lBenomar et al.l 
(l2012al ) to measure the noise characteristics, the 
mode envelope parameters, the mean large separa- 
tion Al/, as well as ep, the phase offset for £ = 
p modes. 

2) After normalizing by the noise background the 
smoothed power spectrum, we identify peaks of 
power for which a null hypothesis is rejected 
at 5% and conside r them as potential modes 
(|Appourchauxll2003[) . 

3) We assigned an identification of £ and n for each 
significant peak that was compatible with Eq|6] (p 
modes) or that follows the expected pattern of 
mixed modes using the asymptotic relation for the 
mixed modes (|Mosser et al.ll20123 ). This gave us a 
first guess for AHi and we checked the rendering in 
a replic ated echelle diagram, as in IBenomar et al.l 
(!2012bh . 

4) A pre-fit (using a MAP approach) to the £ = 1 
peaks provided initial guesses for the matrix A 
(Eql2]). Several iterations were sometimes needed 
between steps 2 and 3, before getting a satisfac- 
tory result. Particular attention was given to the 
number of tt and 7 modes, in order to accurately 
reproduce the avoided crossings. 

5) Guesses on widths of the modes (F) were provided 
by the correlation func tion between Cp and F, found 
bv I White et al.l (|2012| ). Initial heights were defined 
with the help of the envelope h e ight an d Ai/, as 
demonstrated bv lBenomar et all ()2012af ). 

6) The final step was to define priors on the observed 
frequencies. In this paper, we used uniform priors. 
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set after a careful examination of the power spec- 
trum. These priors avoided mode mis-identification 
(in radial order and degree), and sped up the sam- 
pling. 

Steps (3), (4) and (6) rely partly on a visual inspection 
of the power spectra. This inspection is important when 
the spectra show many mixed modes (e.g. giant stars) 
and/or in low signal-to-noise modes. We note however, 
that the asymptotic relations for the p modes and for 
the mixed modes have strong theoretical ground, and 
therefore give confidence to our identifications. More- 
over, simulations (Section 2]) strengthen the validity of 
our method. 

4. TESTS ON SIMULATED SPECTRA 

The method presented above allows us to fit 
the properties of the modes that are statistically 
significant in the power spectrum, namely, 

- The frequencies, heights and widths of the 
observed modes of degree £ = 0,1, 2, 3. 

- The frequencies of the fictitious = 1 7 and 
TT modes as well as their coupling coefficient 



In order to verify the robustness of our method, we 
applied it to spectra generated artificially that mimic the 
main characteristics (modes height, width, frequencies as 
well as the signal-to-noise) of real stellar cases. 

4.1. Frequencies 

We found the observed frequencies (directly observed 
p modes and mixed modes), and the input frequencies 
to agree within 2a with no evidence of biases. It is in- 
teresting to compare the mean large spacing Aiy 
and the mean period spacing AIIi with the input 
values of the simulated spectra because they are 
good indicators of the mass, radius and age of 
the stars. Concerning the inferred C. = 1 7r-mode 
frequencies, the approach works well when the 
density of 7 modes is lower than the density of 
TT modes. For the less complex cases (when the 
number of 7 modes per Ai/ was less than w 1 at 
t'max): the method accurately extracts the depar- 
tures that arise from glitches in the stellar struc- 
ture. However, when the density of the 7 modes 
was higher than that of the tt modes {i.e. during 
the early giant phase, when Aiy w 20 — 35 /iHz) 
the amount of information that comes from the 
power spectrum was not high enough to ensure 
the accuracy. 

Fig. [1] shows the correlation function and the PDFs 
for A J/ and ep, and for AIIi and a-y. This illustrates the 
degree of precision achievable for these global quantities. 
It is clear that the correlations reduce the precision on 
these parameters, and may introduce systematics. The 
observed difference between the true value and the mea- 
sured one is relatively small for ep and Av (within la 
confidence interval), because p modes are directly ob- 
served in the power spectrum. Concerning the £ = 1 7 
modes, our lack of knowledge about and the difficulty 
to infer it (and therefore AIIi) means that the uncertain- 
ties on AIIi and are of about 1 to 5%. In addition. 



we notice slight systematic error, that remains less than 
4%. The eigensolutions of A (EqllJ change slightly as a 
function of the number of implemented modes in the sim- 
ulation. Neglecting the lowest frequency 7 modes 
shifts eigenfrequencies towards lower values and 
overestimates and underestimates AIIi. The ef- 
fect is intrinsic to the adopted methodology and to the 
coupled-oscillators model approach. Ideally, we should 
use the full set of 7 modes, which follow the asymptotic 
relation of the g modes, to ensure a very accurate fit. 
However in practice this is impossible, due to computa- 
tional limitations. In order to limit this bias as much 
as possible, we paid careful attention to this issue by 
checking when eigensolutions are virtually unchanged by 
adding/substracting tt or 7 modes (see Section [3T5|) . Fi- 
nally, we also compared our AIIi and Av values to those 
obtained with asym ptotic relation for mixed modes (cf. 

osser et aI1l2012c[ )) and found that our two approaches 
are in agreement (cf. Fig. [1]). 

With all the assumptions presented above, the 
coupled oscillator model provide a good guide to 
identify and fit the mode frequencies in evolved 
solar-like stars. 

4.2. Height, width and their dependence on the splitting 
and inclination 

The conditions for which the fit is reliable depend on 
whether the assumptions of Section [331 are fulfilled. The 
most important assumption concerns the rotational split- 
ting. Our simulations show that the splitting could not 
be measured in two situations: 

1) When the splitting was smaller than the width of 
the modes and could not be resolved. 

2) When the inclination angle i was below about 30°, 
so that the m ^ components contained only a 
small fraction {< 10%) of the total power of the 
modes and each degree i was well approximated 
by a single Lorentzian. 

When one of these conditions was fulfilled, the mea- 
sured height and width of a given degree has a minimal 
bias. 

To evaluate the bias due to neglecting the rota- 
tional splitting, we fitted a limit-spectrum (with- 
out noise) of an £ = 1 triplet at an inclination of 
90° with a single Lorentzian. This is a worst-case 
scenario because the energy is only on the m = ±1 
components. Fig. [2] shows the results of such a 
fit. The MCMC analysis provides PDF of the fit- 
ted parameters, an d a criterion has to be chosen 
as 'best' fit. Like inlBenomar et al.l ([2009a ^b.) and 
iGruberbauer et al.l ( 2009[ ). we chose the median 
because it is more robust than the mean or the 
maximum of the PDF for noisy data. This sta- 
tistical criterion induces a bias of about 1 % the 
height and the width0- Knowing this, the bias 

^ The maximum of the PDF is a better estimator of 
the height and the width while fitting a limit-spectrum, 
but is not suitable with noisy data because the PDF may 
have several maxima. We also remind that Bayesian ap- 
proaches coupled with MCMC are global optimisation ap- 
proaches and one should rely on the confidence intervals 
rather than on an arbitrary criterion such as the median. 
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Figure 1. Correlation function and the PDFs for the large separation and ep (left), for AIIi and ctj^i (right) measured with the coupled 
oscillators model on a synthetic spectra. Black dashed lines indicate the true values of the parameters. The red dashed line is a linear fit 
of the correlation function. The value retrieved with the asymptotic relation for the mixed modes is shown with a black diamond. The 
correlation introduces quite often biases of a few percent. This bias can be as large as 7% for AEIi. Fixing or ep to their true value 
removes bias. 




Figure 2. Bias on height, width (FWHM, F) and amplitude while fitting a £ = 1 triplet (inclination of 90°) with a single 
Lorentzian as a function of the rotational splitting. The red dotted vertical line indicates the lower limit to visually 
observe the splitting in the power spectrum for F = 1 (/xHz). The bias is more important on height and width for narrow 
modes. Amplitudes are almost not biased. There is a residual bias (ie. when the splitting is of fiHz) due to used 
statistical criterion (median). 



due to neglecting the rotational splitting is negli- 
gible if the rotational splitting remains lower than 
the resolution. When it is not the case, the bias 
is more important for narrow modes, leading to 
an overestimation of the width and an underes- 
timation of the height. However, amplitudes are 
not biased, even when the fit is poor because of 
the anti-correlation between height and width. 

From our simulation, it is clear that our approach is not 
suited for stars whose rotational splitting is much greater 
than the resolution and which have high inclination angle 
(> 30°). 

5. APPLICATION TO FOUR KEPLER STARS 

In this study we have used short-cadence time series 
obtained by the Kepler mission during quarters 5 to 7. 
The corrected flux of each star has been computed fol - 
lowing the procedures described by IGarcia et all (|2011[ ) 
to remove instrumental effects. The power spectral den- 
sity has been obtained using a Lomb-Scargle algorithm 
and normalized following Parseval's theorem (frequency 
resolution of ~ 0.042 /xHz). 



We used three criteria in order to select stars to ana- 
lyze. First, since the approach achieves its full potential 
when the splitting is not resolved or the inclination angle 
is small (cf Section |42|) . we kept only targets for which 
the splitting is not clearly visible in the power spectrum. 
Second, in order to have stars of various evolutionary 
stages and masses we selected targets of different large 
separation {Av from 32 to 65 /zHz) and mean period 
spacings (from 150 to 300 seconds). We kept the four 
stars with the highest signal-to-noise and carried out the 
analysis according to the method described Section 13.51 
This allowed us to measure all the individual properties 
of the modes (frequency, height, width, amplitude). Fig. 
[H sho ws the positions of th ese stars on a AIIi — Av dia- 
gram (iMosser et al.ll2012c| ) which gives us an indication 
on their evolutionary stage and their mass. 

5.1. Frequencies of the modes 

Table [1] contains the main seismic frequency parame- 
ters, computed from the fit. Note that Aiy, ep, AIIi 
and e^^i are determined by fitting a first-order 
polynomial to the frequencies of £ = p modes 
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Figure 3. AIIi — Ai' diagram for the stars of tlie present study 



(black crosses) superimposed on tiie results from iMosser et al.l 
Il2012d ) for red giant stars (gray diamonds). Tlie dotted line in- 
dicates wlien Nt^ = Nj. Solid colored lines are models for masses 
between I.OMq and I.SMq. 

and £ = 17 modes, at each iterative step of the 
MCMC proces^. Glitches inside stars introduce 
oscillations on Aiy and AIIi that can slightly dis- 
tort their determination, depending on the ob- 
served mode range. Therefore, to compare our 
values with theoretical Ai/ and AIIi one needs to 
use the same range of frequencies, as specified in 
the Tables [5]l8l These tables also contain the in- 
dividual mode parameters of the analy zed stars. 

Scaling relations, were already discus sed by lEddingtonl 
([1913) for Cepheids (pBelkac em |2012| ) . Howeve r scal- 
ing re lations have been in t roduc ed in the 80's by lUlrichI 
()1986[) and iBrown et al.l (I1991D before being u sed for 
solar-like oscillators (jKieldsen fc Bedding! |1995() . We 
used them to determine the masses and radii in Table 
[21 Note that scaling relations may suffer from system- 
atics in some part of the HR diagram (such as for red 
clump star) becau se they are calibrated using the Sun 
(|Miglio et al.ll2012D. Effe c tive t empe ratures Teg in Table 

ire from iBruntt et al.l ()2012D and iPinsonneault et all 

)12). 

Figs. [4]l7] show echelle diagrams (a) and repli- 
cated echelle diagrams (b) for the four stars. 
These diagrams show all fitted frequencies (p, 
TT, 7 modes, and mixed modes). KIC 6442183 
(Fig. H) and KIC 11026764 (Fig. [5]) have only 
a few avoided crossings and are therefore early 
subgiants. These avoided crossings can be seen 
as departures from the asymptotic relation for 
the p modes in the echelle diagram. 

Interestingly, in addition to two clear £ — 1 avoided 
crossings, KIC 1102 6764 (also known as Gemma e.g. 
iMetcalfe et al.l |2010[ ). shows one clear £ = 2 avoided 
crossing (Fig. [H), which was expected from models of 
this star but had not been observed before. Their po- 
tential interest is discussed Section 15.21 The coupled 
oscillator model extracts very precisely the position of 
the TT modes allowing an accurate determination of the 

^ We therefore obtain the PDF for these asymptotic 
parameters from which we derive the median and the un- 
certainty. 



frequency separation of each mixed modes from the tt 
modes home ridge, denoted d in the following. This per- 
mits us to quantify the effect of the degree of mixture, 
on the individual modes parameters, as discussed below 
in Section l573l 

KIC 12508433 (Fig. ^ and KIC 11771760 (Fig. ^ 
are very interesting because contrary to KIC 6442183 
and KIC 11026764, the density of 7 modes exceeds the 
density of tt modes, showing that these stars are more 
evolved. Their positions in the AIIi-Ai/ diagram indicate 
that they are at the transition be tween the subgiant and 
giant phases (jMosser et al.ll2013l in prep.). 

Several £ = 3 modes are detected in KIC 6442183 
and KIC 12508433. These modes behave as pure p 
modes because the evanescent zone is larger. Ex- 
amination of the echelle diagram of KIC 11026764 with 
a strong smoothing, shows a slight power close to the tt- 
mode ridge (Fig. [S|). These possible £ = 3 modes were 
not fitted because they are below the detection thres hold 
with which we defined the initial guesses (cf. Section [331 
step 2). It is likely that more observations will allow an 
accurate fit to these i = S modes, providing additional 
constraints to models. 

5.2. Testing the asymptotic relation for gravity modes 
using the 7 modes 

Measuring gravity-modes of different degrees can 
greatly improve the constraints on the properties of stel- 
lar cores. From EqlTl we can easily show that An2 = 
A^l/^/3, which means that the density of ^ = 2 modes 
is expected to be greater than ior £ = 1. More interest- 
ingly, we also obtain. 



£7,2 ~ £7,1 



1/2 (radiative core) 
(convective core) 



(10) 



Therefore, if both phase offsets for £ = 1 and ^ = 27 
modes could be measured, their difference can only take 
two values, independent from the actual phase offs et a-y . 
This removes the degeneracy on (see Section 13.4. ip 
and potentially provides a direct insight on the nature of 
the core, independent to the stellar models. 

Subgiant and RGB stars should not have a convective 
core as Hydrogen is depleted in the core. Therefore con- 
vection is not required to transport energy and we can 
reasonably assume that e-y,2 — e-y,i — 1/2. However, Eq[7l 
is a first-order approximation and assumes that the radial 
order of the gravity modes is high. In early subgiants, 
the observed £ = I avoided crossings correspond to low 
order gravity modes. This implies that AH2 may not 
scale with AHi/VS and that aj,e = + £/2 may no 
longer be valid. Therefore, measuring e-y.2 — e-y.i allows 
us to probe the accuracy of the asymptotic relation for 
the gravity modes. 

As discussed in Section 13.4. 2[ the mode bumping of 
£ = 2 mixed modes can hardly be observed as the cou- 
pling between the pressure and gravity modes is weaker 
than for £ = 1. Nonetheless, among the four analyzed 
stars, KIC 11026764 and KIC 12508433 show clear £ = 2 
bumped modes. In order to quantify approximately e-y.2i 
we performed a fit of the measured £ = 2 frequencies 
(output from the MC MC fit) in a similar fashion to 
iBenomar et al.l (l2012b( ) and derived the expected fre- 
quencies of £ = 2 7-mode ^■y,2, causing the most sig- 
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Figure 4. EchcUc diagrams of the fitted frequencies (a) and replicated echelle diagrams (b) for KIC 6442183. Orange, gray, yellow 
colored dots represent £ = 0,2,3 pure p modes. £ = 1 fictitious tt modes are blue while £ = I observed mixed modes are 
red. Horizontal dashed-red lines show where i = 1 avoided crossing occur (7- mode frequencies). All uncertainties are shown at 3a. 




Figure 6. Echelle diagrams of the fitted frequencies (a) and replicated echelle diagrams (b) for KIC 12508433. Same legend as Fig. |3] 
The circle dot indicates the expected frequency of an £ = 1 mode but not statistically significant. 
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Figure 7. EchcUe diagrams of tiic fitted frequencies (a) and replicated eclielle diagrams (b) for KIC 11771760. Same legend as Fig. |4] 
The sharp feature on the £ = 1 it modes is presumably an artifact. The circle dot indicates the expected frequency of an f = 1 mode but 
not statistically significant. 
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Figure 8. KIC 11026764 spectrum showing the i = 2 mixed 
modes and the £ = 2 7-mode {^■^1^2)- Closest 1 = and 1=1 
modes are also represented. Thick-colored horizontal lines repre- 
sent confidence intervals at 2cf. 

nificant mode bumping. 

For KIC 11026764, we obtained v^^2 = 768.70 ± 0.80 
/iHz. Using EqlT] and the measured AIIi this leads to 
e^,2 = 0.49 ±0.20. We have e^_i = 0.59 ±0.20 (see Table 
[TJ whieh gives a differenee e.'y,i=2 ~ £7,^=1 = 0-10 ± 0.20. 
This value is rather low and suggest that v^e are 
not in the asymptotic regim^^. 

For KIC 12508433, whieh is more evolved, the main 
mode bumping is at i'^_2 = 737.74 ± 4.50 /iHz and we 
obtain e^^2 = 0.91 ± 0.02 which leads to £7,^=2 — £7.^=1 = 
0.40 ± 0.06. This is more consistent with the asymptotic 
relation for the gravity modes than for KIC 11026764. 

5.3. Heights, widths and amplitudes of the modes 

While pulsation frequencies provide constraints on the 
mass, radius and evolutionary stage of stars, the widths, 
heights and amplitudes of the modes may carry informa- 
tion on the mode inertia and on the non-adiabatic 
effects such as the radiative damping. Measuring 
them is therefore essential to constrain stellar models and 
to better understand the role of convection in the exci- 
tation mechanism of the modes. 

Due to the long observations achieved with Kepler^ 
many mixed modes are well resolved and an accurate 
fit is possible. Modes a posteriori not statistically signif- 
icant are discarded in the following analysis. 

It is well known that width and height of modes 
are an ti-correlated (e . e;. [ Appourchaux et al.l 
119981 : iBenomar et al.l l2009an . This anti- 

correlation greatly reduces the precision of their 
measure, whereas their product, which defines 
the mode amplitude, is much better constrained. 
The relationship between height, width and am- 
plitude is well defined, therefore in the follow- 
ing we will focus mainly on width and amplitude, 
without loss of information. 

Fig. [9l and [TOl show the Full Width at Half Max- 
imum (FWHM) Ti and the amplitudes Ai (with 
A( ~ ^/-kHiTi). In each case, we show the actual 

* Note that n-, 1 = 4, which also suggests that we are not 
in the asymptotic regime. 



value as a function of frequency (left column) and 
a normalized value as a function of the distance 
d, to the £ = 1 TT-home ridge. The normalized val- 
ues were calculated using the local, linearly interpolated 
I = {) values as a reference. This allows us to compare 
the relative width (in /^Hz) and energy (ppm) of 
£ =1 modes for all the stars. 

We note that the width for € = 1 is generally lower than 
for £ = and drops at each avoided crossing. These drops 
are mainly an effect of the coupling between tt and 7 
modes. The more gravity-like modes are narrower which 
can also be seen clearly in plots of FWHM versus d/ /S.v. 
A similar effect is seen more weakly in the heights 
of modes. However, when the mode FWHM ap- 
proaches the spectrum resolution the height un- 
certainty becomes great and the stochastic exci- 
tation of the mode prevents us from providing 
a reliable measure. The effect of the coupling is 
clearest for KIC 11771760. The number of mixed 
modes is greater. In the T{d//S.v) and A{d//S.v) 
diagrams, the normalized quantities follow a bell 
shape, centered on approximately showing that 
modes with mainly a p-mode character are the 
more intense, as expected. 

In addition to the mixing, geometrical effects and limb 
darkening modify the relative values of heights and am- 
plitudes, and therefore the visibilities of each degree. 
Before quantitative comparison with models, one would 
need to disentangle these effects. If we assume that vis- 
ibilities of all € = 1 modes are identical, we can esti- 
mate the visibility by measuring the relative heights or 
amplitudes at the home ridge, meaning when d « 0, 
where mixed modes behave as p modes. To do this, 
we fitted the measured height and amplitudes with a 
parabolic function. However, a stable fit was not al- 
ways achievable, especially for the early subgiants. In 
this case, we fitted a constant around the maximum of 
the function d/Ai/, by including only values that satisfies 
|(i/Aj/| < 0.05. Table |3] and |4] summarize the obtained 
visibilities for all non-radial modes in height and ampli- 
tude. For conipariso n purpose, solar height visibilities 
(iSalabert et al.ll2011f) and for average RGB stars visibil- 
ities (jMosser et al.ll2012bl) are also shown. 

The £ = 1 relative heights suffer from high uncertain- 
ties. However the early subgiant (KIC 6442183 and KIC 
11026764) visibilities are rather high compared to the 
Sun and RGB stars while for stars in the late subgiant 
phase (KIC 12508433 and KIC 11771760), values become 
comparable to RGB stars. It is clear that £ = 2 heights 
relative to £ = 0, are substantially higher than in the 
Sun or RGB for KIC 11026764 and KIC 12508433. This 
has been already noticed in several su bgian t stars (e.g. 
iDeheuvels et al.l (|2010f ): iMathur et al.l (|2013| )). although 
with greater uncertainties. The £ = 3 visibilities are 
consistent with solar visibilities for both subgiants and 
RGBs. According to iBallot et all ()20lH) the visibility 
factors depend on the effective temperature and on the 
metallicity. However, these cannot entirely explain the 
discrepancy between main-sequence, subgiant and giant 
stars. 

6. DISCUSSION 

Precise measures of all the properties of the os- 
cillations in early subgiants (i.e. when <C A^tt) 
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Figure 9. £ = (orange) and £ = 1 (blue) widths (FWHM, F). Vertical red dashed lines indicate 7 modes positions. Vertical black dotted 
lines indicate tt modes positions. Red data points arc not statistically significant in the power spectrum. When the i = I mode is far from 
the TT mode frequencies, width drops. 
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Table 1 

Tabic of computed seismic parameters for the four stars. 



KIC I'max (mHz) IS.V (/jHz) AIIi (sec) e^.i (5i/oi (/iHz) hv^2 (mHz) Si^is (aiHz) 

6442183 1160 ±4 65.07 ± 0.09 1.467 ±0.021 325 ± 18 0.05 ±0.11 4.2 ±0.2 5.60 ±0.20 7.00 ± 0.23 

11026764 880 ±5 50.17 ±0.09 1.453 ±0.025 301 ± 16 0.59 ± 0.20 1.45 ± 0.6 4.14 ±0.08 / 

12508433 773 ± 5 45.18 ± 0.05 1.407 ±0.008 157.5 ± 0.8 0.524 ± 0.034 -0.7 ± 0.2 4.02 ±0.10 10.91 ±0.10 

11771760 522 ±7 32.34 ± 0.09 1.206 ±0.036 177 ± 5 0.173{+0.30 - 0.24) 1.28 ± 0.38 3.00 ± 0.09 / 
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Figure 10. £ = (orange) and £ = I (blue) Amplitude. Vertical dashed lines indicate 7 modes positions. Vertical black dotted lines 
indicate n modes positions. Red data points are not statistically significant in the power spectrum. When the £ = I mode is far from the 
TT mode frequencies, amplitude drops. On right, the orange dashed line is a polynomial fit performed to measure the visibilities. 
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Table 2 

Effective temperature and derived seismic radius and mass of tlic 
four stars using scaling relations. Typical uncertainties on mass 
and radius are a few percent. 



KIC 


Teff (K) 


M (Mq) R (Rq) 


6442183 
11026764 
12508433 
11771760 


5740 ± 60 
5722 ± 57 
5257 ± 64 
6075 ± 70 


~ 0.94 ~ 1.60 
~ 1.14 ~2.02 
~ 1.03 ~2.10 
~ 1.51 ~2.98 


Table 3 

Height visibility factors for £ = 1,2,3 modes. 


KIC 


He = 1/Hi = 


He = 2/ He = He = 3/ He = 


6442183 
11026764 
12508433 
11771760 


1.88 ±0.15 
2.04 ±0.21 
1.50 ±0.31 
1.26 ±0.10 


0.584 ±0.055 0.081 ±0.015 
0.773 ± 0.081 / 
0.581 ± 0.050 0.082 ± 0.016 
0.776 ± 0.088 / 


Sun 
RGB 


1.53 
1.35 


0.56 0.034 
0.64 0.071 


Table 4 

Amplitude visibility factors for £ = 1,2,3 modes. 


KIC 


= 1/Ae = 


Ae = 2/Ae = Ae = 3/Ae = 


6442183 
11026764 
12508433 
11771760 


1.27 ±0.15 
1.37 ±0.12 
1.18 ±0.09 
1.09 ±0.14 


0.751 ±0.061 0.285 ±0.032 
0.875 ± 0.082 / 
0.768 ±0.063 0.272 ±0.034 
0.809 ± 0.075 / 



are achieved with the coupled oscillator model 
without any problem. Inacurate reconstruction 
of the TT-mode and 7-mode spectrum (artifacts) 
may yet happen in late giants with the cur- 
rent methodology. However, because the more 
evolved stars have a very dense 7-mode spectra 
{N^ ^ iV^r) a rich mixed mode spectrum is ob- 
served that might provide more stringent con- 
straints on the internal structure for giants than 
for subgiants. Further work might therefore focus 
on late giants if one wants to, for example, apply 
inversion methods to constrain the internal stel- 
lar structure using the widths and amplitudes of 
mixed modes. 

In this study, we restricted our analysis to the simplest 
subgiants, without visible rotation, with a low gravity- 
mode density and at high signal to noise. However, in 
many cases, observed subgiants and giants stars present 
a clear rotational signature and often have a radial differ- 
ential rotation, implying important variations of splitting 
from one ridge to another. Modes width could then be 
overestimated in the case of a non-negligible splitting. 
With stellar cores, ro tating faster than the surface {e.g. 
iDeheuvels et al.ll2012[ ). the bias will be more significant 
for modes far from the p modes home ridge. 

A more sophisticated analysis, which has still to be de- 
veloped, would be necessary in order to measure precisely 
height, widths and amplitudes, but also rotational (and 
possible magnetic) splittings of all modes. At the mo- 
ment, global approaches that use the coupled oscillator 
model require MCMC fits to provide reliable results and 



are complex. This study is at the edge of what currently 
could be possibly done. The main obstacle is the compu- 
tational time, which increases quickly with the number 
of degrees of freedom of the fit. To overcome this diffi- 
culty, a local fitting approach, focusing on the stars with 
long-continuous time-series and with the highest signal- 
to-noise is probably more suited. 

In this paper, we also quantified for the first time the 
effect of the coupling between pressure and gravity modes 
on widths, heights and am plitudes. Q ualitatively our re- 
sults are in agreement with lDupret et a l. (2009;), however 
a quantitative comparison would certainly give a new in- 
sight on stellar interiors. 

Interestingly KIC 11026764 and KIC 12508433 have 
clear bumped £ = 2 modes. As for £ = 1 modes, this 
mode bumping is a signature of their mixed character. 
However, £ = 2 mixed modes have a weaker coupling 
than in ^ = 1 modes and therefore avoided crossings are 
not easy to detect in the echelle diagram. For KIC 
11026764, the two mixed modes are measurable. 
But for KIC 12058433 only a bumping on one 
of them is visible, the second mixed mode being 
either too faint or not resolved. 

A combined measure of = 1 and i = 2 mixed modes 
allows us to better constrain stellar core properties. Here, 
we showed that the phase offset of the gravity modes can 
help in the asymptotic regime to constrain the nature 
of the core (convective or radiative) . However this phase 
offset is very sensitive to the inaccuracy of the asymptotic 
relation for the gravity modes. When the nature of the 
core is unambiguous, it can be used to delimit the validity 
of the asymptotic theory for the gravity modes. In giants 
and subgiants, no fusion reactions occur in the core and 
therefore it is necessarily radiative. 

In a subgiant such as KIC 11026764, the asymptotic 
relation for the gravity modes is a poor approximation 
because it is a first-order approximation and is valid only 
for high radial order. In RGBs, the mode order at max- 
imum of amplitude is rather high, such that the asymp- 
totic relation for the gravity modes is fairly accurate. 
KIC 12058433, being an early red giant, has a phase off- 
set consistent with theoretical expectations. 

7. CONCLUSION 

Mixed modes are difficult to identify and fit because 
their frequencies do not follow the asymptotic relations 
for pressure or gravity modes. However, as shown in 
the present paper, they are not randomly distributed in 
frequency but rather follow a distinctive pattern. This 
pattern is a sum of several avoided crossings, each being 
due in subgiants to the coupling of pressure modes with 
a single gravity mode. 



Methods and tools to i dentify mixed modes have 



been already proposed (e.g. 


Beddinsll2012l:IMosser et al.l 


l2012cl:lBenomar et al.ll2012b 


). But among them the cou- 



pled harmonic oscillator model is the only one that re- 
lies on individual pressure and gravity mode frequencies 
(rather than global, average quantities such as Aiy or 
AHf) and this allows to fit precisely the observed power 
spectrum. However, with this model a Bayesian ap- 
proach is necessary because it is very complex and has 
too many degrees of freedom. Without guidance the 
model could easily lead to non-physical results. 
In the present work, we used a set of priors that de- 
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Table 5 

Observed mode parameters and their uncertainties for KIC 6442183. 



KIC 6442183 - Mode range: [700 - 1500] ^Hz 



i 


7 TT \ 


7 \ 


7 — \ 

err- (u) 


TT? TT \ 

r (/iHz) 


err- (1 ) 


err+(l ) 


K 7 \ 

A (ppm) 


err- (A) 







743.47 


0.11 


0.11 


0.64 


0.14 


0.13 


2.59 


0.48 


0.46 







0.23 


0.16 


0.88 


0.22 


0.24 


2.51 


0.45 


0.45 





873.85 


0.10 


0.10 


0.83 


0.14 


0.17 


4.06 


0.39 


0.39 





937.14 


0.10 


0.10 


0.82 


0.11 


0.14 


4.73 


0.37 


0.39 





1001.04 


0.09 


0.09 


0.74 


0.09 


0.11 


6.06 


0.49 


0.52 





1 065 63 


0.09 


0.09 


0.79 


0.10 


0.11 


6.46 


0.42 


0.45 





1130.49 


0.06 


0.06 


0.69 


0.08 


0.09 


9.61 


0.53 


0.55 





1195.41 


0.06 


0.06 


0.66 


0.08 


0.08 


9.87 


0.55 


0.60 





1 260 23 


0.08 


0.09 


1.13 


0.13 


0.14 


7.85 


0.41 


0.43 





1 325 63 


0.14 


0.13 


1.46 


0.19 


0.21 


6.58 


0.39 


0.40 





1391.92 


0.22 


0.22 


2.07 


0.28 


0.37 


4.76 


0.33 


0.33 





1458.12 


0.53 


0.47 


3.37 


0.49 


0.58 


3.31 


0.29 


0.29 


\ 


745,74 


5.61 


0.24 


0.25 


0.16 


0.29 


1.74 


0.62 


0.62 


1 


783.04 


0.30 


0.31 


0.64 


0.22 


0.20 


3.18 


0.54 


0.53 


1 


840.74 


0.09 


0.09 


0.62 


0.18 


0.22 


4.26 


0.49 


0.49 


1 


901.28 


0.09 


0.10 


0.70 


0.16 


0.17 


4.73 


0.45 


0.48 


1 


960.29 


0.10 


0.10 


0.82 


0.15 


0.15 


5.69 


0.48 


0.53 


1 


1002 93 

_L \J\J t\jKJ 


0.06 


0.05 


0.26 


0.09 


0.18 


6.17 


0.80 


0.89 


\ 


1037 59 


0.07 


0.07 


0.64 


0.13 


0.13 


8.07 


0.60 


0.67 


\ 


1 096 78 


0.05 


0.05 


0.51 


0.10 


0.12 


9.98 


0.78 


0.89 


1 


1159.92 


0.05 


0.05 


0.55 


0.10 


0.10 


13.15 


0.90 


1.17 


1 


1224.10 


0.07 


0.06 


0.85 


0.12 


0.12 


11.32 


0.68 


0.73 


1 


1288.63 


0.09 


0.09 


1.23 


0.17 


0.16 


8.99 


0.51 


0.55 


1 


1 353 98 


0.12 


0.12 


1.71 


0.22 


0.23 


7.78 


0.40 


0.45 


1 


1419.44 


0.14 


0.13 


1.67 


0.30 


0.34 


5.82 


0.37 


0.39 


2 


737.14 


2.54 


1.11 


0.58 


0.13 


0.15 


1.92 


0.39 


0.40 


2 


801.89 


0.32 


0.70 


0.90 


0.25 


0.28 


2.01 


0.36 


0.35 


2 


867.84 


0.20 


0.19 


0.78 


0.18 


0.22 


2.91 


0.33 


0.33 


2 


931.41 


0.14 


0.15 


0.82 


0.15 


0.18 


3.58 


0.31 


0.34 


2 


996.73 


0.10 


0.10 


0.52 


0.10 


0.12 


3.79 


0.34 


0.35 


2 


1059.92 


0.08 


0.07 


0.71 


0.11 


0.13 


4.71 


0.32 


0.33 


2 


1124.66 


0.06 


0.06 


0.63 


0.08 


0.09 


6.82 


0.40 


0.42 


2 


1189.77 


0.07 


0.07 


0.72 


0.09 


0.11 


7.88 


0.41 


0.44 


2 


1254.54 


0.11 


0.11 


1.20 


0.15 


0.16 


6.67 


0.36 


0.35 


2 


1321.35 


0.14 


0.14 


1.25 


0.21 


0.24 


4.79 


0.33 


0.33 


2 


1386.89 


0.21 


0.20 


2.30 


0.36 


0.43 


4.08 


0.26 


0.28 


2 


1453.09 


0.38 


0.38 


3.08 


0.55 


0.65 


2.63 


0.22 


0.22 


3 


1085.64 


1.01 


0.32 


0.76 


0.08 


0.09 


2.19 


0.20 


0.22 


3 


1150.25 


0.17 


0.17 


0.76 


0.11 


0.12 


2.92 


0.26 


0.28 


3 


1215.67 


0.18 


0.15 


0.66 


0.15 


0.16 


2.51 


0.28 


0.29 


3 


1280.75 


0.39 


0.35 


1.42 


0.28 


0.31 


2.34 


0.27 


0.27 


3 


1347.01 


0.37 


0.41 


1.41 


0.32 


0.37 


1.65 


0.22 


0.22 


3 


1413.29 


0.77 


0.71 


2.50 


0.36 


0.43 


1.31 


0.15 


0.16 



scribe as best as possible our current knowledge, achiev- 
ing proper balance between prior and likelihood. The 
coupled oscillator model has been embedded in a MCMC 
fitting algorithm in order to both identify and fit the in- 
dividual gravity, pressure and t ~ 1 mixed modes of sub- 
giant stars. Thanks to this, we were able to fit precisely 
all the properties (width, height and frequency) of the 
observed modes for four subgiants of different mass and 
evolutionary stage. The coupled oscillator model also 
permits us to infer the frequencies of the gravity-modes 
and therefore to measure accurately the period spacing 
of the gravity-modes. 

The coupled oscillator model allows us to infer 
where the 1=1 pressure modes would be if pres- 
sure and gravity modes were uncoupled (referred 
to as the home ridge). These mode frequencies in 
return, permit us to measure the actual distance 
between the home ridge and the observed mixed 
modes, which is a proxy of the degree of mix- 
ing. We therefore measured the effect of the mix- 
ing on the main characteristics of the modes and 
found that the more gravity-like the modes, the 



lower the widths and amplitud es. This qualita- 
tively confirms the results of Dupret et al. (2009). 
We stress that the stochastic nature of the ex- 
citation of the solar-like oscillations complicates 
measurements of mode heights. However, analy- 
sis of longer time-series should allow us better to 
constrain this parameter thanks to an increase of 
the frequency resolution. 

Finally, in addition to ^ = 1 mixed modes, we also 
clearly identified an £ = 2 mixed mode in KIC 11026764 
and showed that we are not in the asymptotic regime for 
gravity modes. 
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Table 6 

Observed mode parrameters and their uncertainties for KIC 11026764. 











KIC 11026764 - 


Mode range: [550 - 


- 1100] /^Hz 








£ 






G.TV {^l* 


) ^ v/^J^^; 


prr (V\ 
L.I 1 -yi ) 






GT'T' {-^^ 


e; f + 





574.68 


0.11 


0.16 


0.39 


0.14 


0.23 


2.98 


0.47 


0.52 





625.01 


0.10 


0.11 


0.40 


0.10 


0.13 


3.92 


0.43 


0.47 





674.53 


0.09 


0.10 


0.52 


0.14 


0.24 


3.76 


0.44 


0.47 





723.48 


0.15 


0.14 


0.99 


0.20 


0.25 


5.11 


0.50 


0.49 





773.00 


0.11 


0.09 


0.85 


0.16 


0.19 


4.59 


0.37 


0.42 





823.21 


0.07 


0.07 


0.76 


0.11 


0.13 


7.93 


0.51 


0.52 





873.75 


0.07 


0.07 


0.78 


0.10 


0.12 


8.64 


0.53 


0.56 





924.06 


0.10 


0.10 


0.90 


0.14 


0.16 


8.27 


0.57 


0.55 





974.56 


0.15 


0.14 


1.22 


0.17 


0.21 


6.99 


0.46 


0.48 





1025.49 


0.17 


0.17 


1.69 


0.26 


0.34 


5.61 


0.40 


0.41 





1076.10 


0.52 


0.42 


2.92 


0.55 


0.75 


3.66 


0.36 


0.37 


1 


694.12 


0.06 


0.06 


0.37 


0.13 


0.16 


5.14 


0.58 


0.72 


1 


723.91 


0.33 


2.47 


0.03 


0.02 


0.16 


2.15 


1.16 


1.44 


1 


755.10 


0.10 


0.11 


0.83 


0.17 


0.16 


6.45 


0.53 


0.63 


1 


799.62 


0.08 


0.08 


0.78 


0.14 


0.14 


8.91 


0.62 


0.73 


1 


847.39 


0.07 


0.07 


0.73 


0.17 


0.14 


9.85 


0.74 


1.10 


1 


894.03 


0.06 


0.06 


0.67 


0.13 


0.14 


11.56 


0.77 


0.92 


1 


925.83 


0.06 


0.06 


0.39 


0.12 


0.12 


7.99 


0.81 


0.98 


1 


954.09 


0.08 


0.08 


0.96 


0.14 


0.15 


10.67 


0.64 


0.70 


1 


1000.41 


0.08 


0.08 


1.02 


0.18 


0.21 


8.96 


0.60 


0.67 


1 


1050.11 


0.14 


0.14 


1.70 


0.33 


0.37 


6.20 


0.41 


0.46 


2 


570.32 


0.18 


0.14 


0.37 


0.15 


0.23 


2.44 


0.45 


0.50 


2 


620.64 


0.11 


0.09 


0.42 


0.12 


0.15 


3.48 


0.40 


0.43 


2 


669.49 


0.15 


0.13 


0.41 


0.15 


0.23 


3.03 


0.41 


0.41 


2 


718.92 


0.22 


0.22 


1.10 


0.24 


0.30 


4.76 


0.42 


0.42 


2 


765.85 


0.14 


0.21 


0.36 


0.16 


0.21 


2.67 


0.65 


0.59 


2 


771.41 


0.26 


0.27 


0.83 


0.18 


0.24 


4.04 


0.36 


0.38 


2 


818.67 


0.10 


0.10 


0.79 


0.13 


0.16 


6.93 


0.43 


0.47 


2 


868.89 


0.09 


0.10 


0.87 


0.12 


0.13 


7.97 


0.45 


0.49 


2 


919.89 


0.10 


0.10 


1.05 


0.15 


0.18 


7.98 


0.48 


0.49 


2 


970.15 


0.10 


0.11 


1.16 


0.17 


0.20 


6.26 


0.38 


0.42 


2 


1020.41 


0.18 


0.20 


1.50 


0.30 


0.32 


4.89 


0.36 


0.38 


2 


1072.09 


0.45 


0.50 


2.74 


0.56 


0.69 


3.48 


0.32 


0.32 
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Table 7 

Observed mode parrametcrs and their uncertainties for KIC 12508433. The mode with an asterisk is non-statistically significant. 



KIC 12508433 - Mode range: [590 - 1050] fiHz 





u (fiUz) 


err^(u} 


err- (u) 




err- (1') 


err_|_(l') 


A (ppm) 


err -(A) 


err_|_(A) 





606.10 


0.05 


0.06 


0.41 


0.09 


0.12 


4.53 


0.45 


0.46 





649.96 


0.06 


0.06 


0.47 


0.08 


0.10 


6.31 


0.56 


0.57 





694.41 


0.05 


0.05 


0.45 


0.05 


0.06 


7.89 


0.61 


0.64 





739.42 


0.05 


0.05 


0.44 


0.05 


0.06 


11.19 


0.82 


0.90 





784.16 


0.05 


0.05 


0.42 


0.06 


0.07 


13.00 


0.83 


0.91 





829.08 


0.06 


0.06 


0.54 


0.07 


0.08 


10.42 


0.68 


0.72 





874.25 


0.07 


0.07 


0.69 


0.08 


0.10 


9.20 


0.60 


0.65 





919.82 


0.15 


0.15 


1.23 


0.16 


0.17 


5.68 


0.45 


0.49 





965.75 


0.21 


0.23 


1.58 


0.38 


0.44 


4.02 


0.45 


0.45 





1012.57 


0.22 


0.13 


1.97 


0.51 


0.63 


2.91 


0.46 


0.45 


1 


627.14 


0.06 


0.06 


0.38 


0.08 


0.09 


5.31 


0.58 


0.61 


1 


657.25 


0.07 


0.07 


0.26 


0.05 


0.08 


7.51 


0.78 


0.72 


1 


679.15 


0.05 


0.06 


0.37 


0.09 


0.09 


7.01 


0.72 


0.78 


1 


713.54 


0.05 


0.05 


0.35 


0.07 


0.10 


8.84 


0.83 


0.89 


1 


740.66 


0.07 


0.06 


0.25 


0.11 


0.11 


8.42 


1.20 


1.49 


1 


768.29 


0.04 


0.04 


0.32 


0.04 


0.06 


13.98 


0.97 


0.90 


1 


806.14 


0.04 


0.04 


0.32 


0.04 


0.05 


15.21 


0.92 


0.92 


1 


836.05 


0.04 


0.04 


0.31 


0.05 


0.07 


11.46 


0.98 


0.92 


1 


860.47 


0.06 


0.05 


0.51 


0.09 


0.12 


9.04 


0.77 


0.81 


1 


899.99 


0.06 


0.07 


0.59 


0.15 


0.14 


9.14 


0.80 


1.11 


1 


942.31 


0.14 


0.14 


1.27 


0.25 


0.31 


6.02 


0.54 


0.55 


1* 


974.11 


0.14 


0.15 


0.01 


0.01 


0.01 


0.86 


0.41 


0.47 


1 


1000.05 


0.31 


0.33 


1.70 


0.57 


0.64 


4.37 


0.54 


0.54 


2 


602.73 


0.18 


0.19 


0.40 


0.10 


0.14 


3.30 


0.42 


0.46 


2 


646.12 


0.14 


0.12 


0.49 


0.10 


0.12 


4.86 


0.46 


0.47 


2 


690.61 


0.14 


0.14 


0.44 


0.08 


0.09 


5.88 


0.53 


0.56 


2 


735.95 


0.06 


0.06 


0.47 


0.06 


0.07 


8.65 


0.57 


0.58 


2 


779.60 


0.05 


0.05 


0.43 


0.06 


0.08 


9.90 


0.62 


0.64 


2 


825.20 


0.07 


0.08 


0.58 


0.08 


0.09 


8.64 


0.53 


0.57 


2 


870.92 


0.08 


0.09 


0.68 


0.10 


0.11 


7.12 


0.43 


0.48 


2 


915.61 


0.18 


0.19 


1.34 


0.21 


0.24 


5.27 


0.40 


0.42 


2 


961.35 


0.28 


0.56 


1.40 


0.44 


0.50 


3.08 


0.35 


0.35 


2 


1006.98 


0.55 


0.57 


1.92 


0.52 


0.70 


2.39 


0.34 


0.34 


3 


752.04 


0.23 


0.27 


0.44 


0.05 


0.06 


3.23 


0.37 


0.39 


3 


797.26 


0.13 


0.16 


0.44 


0.08 


0.09 


3.33 


0.38 


0.40 


3 


842.53 


0.18 


0.24 


0.54 


0.10 


0.11 


2.66 


0.33 


0.36 


3 


887.44 


0.30 


0.38 


0.90 


0.13 


0.13 


2.50 


0.30 


0.31 
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Table 8 

Observed mode parameters and their uncertainties for KIC 11771760. The mode with an asterisk is non-statistically significant. 



KIC 11771760 - Mode range: [350 - 700] ^Hz 







err^(u) 


erT_ [u) 


r (/iHz) 


err_(l') 


err-|_(l') 


A (ppm) 


erT_(A) 


err_|_(A) 





361.26 


0.08 


0.12 


1.02 


0.23 


0.24 


8.25 


1.07 


1.06 





394.35 


0.18 


0.18 


1.13 


0.24 


0.31 


5.49 


0.66 


0.58 





426.66 


0.10 


0.11 


0.81 


0.10 


0.11 


10.32 


0.67 


0.62 





457.27 


0.11 


0.10 


1.12 


0.11 


0.11 


14.38 


0.68 


0.80 





488.96 


0.17 


0.16 


1.01 


0.17 


0.23 


8.39 


0.75 


1.08 





521.59 


0.13 


0.16 


1.06 


0.10 


0.12 


12.20 


0.61 


0.55 





554.69 


0.10 


0.08 


1.05 


0.11 


0.13 


14.95 


0.70 


0.79 





586.88 


0.18 


0.40 


2.19 


0.26 


0.18 


7.78 


0.67 


0.79 





619.55 


0.28 


0.17 


1.45 


0.18 


0.18 


10.38 


0.85 


0.71 





651.92 


0.11 


0.47 


1.20 


0.13 


0.13 


8.39 


0.54 


0.82 


1 


364.83 


0.24 


0.21 


1.11 


0.24 


0.37 


8.22 


0.68 


0.95 


1 


374.66 


0.27 


0.21 


0.87 


0.41 


0.45 


7.31 


1.32 


1.09 


1 


393.19 


0.50 


0.50 


0.37 


0.25 


0.29 


1.58 


0.68 


0.79 


1 


405.82 


0.12 


0.11 


0.74 


0.29 


0.26 


6.36 


0.88 


0.71 


1 


424.91 


0.20 


0.14 


0.40 


0.28 


0.20 


1.76 


1.08 


1.57 


1 


439.69 


0.14 


0.14 


0.82 


0.17 


0.17 


9.17 


0.95 


0.88 


1 


459.56 


0.21 


0.16 


0.64 


0.52 


0.20 


2.04 


1.00 


0.98 


1 


474.92 


0.10 


0.09 


0.73 


0.14 


0.15 


10.65 


0.93 


1.13 


1 


495.47 


0.08 


0.09 


0.55 


0.13 


0.13 


8.14 


0.93 


0.82 


1 


512.92 


0.10 


0.10 


0.88 


0.19 


0.18 


12.14 


1.03 


1.22 


1 


535.97 


0.12 


0.11 


1.07 


0.12 


0.10 


13.13 


0.69 


0.68 


1* 


557.02 


0.03 


0.05 


0.20 


0.14 


0.20 


0.80 


0.47 


0.75 


1 


575.10 


0.20 


0.18 


1.92 


0.14 


0.17 


13.47 


1.14 


1.37 


1 


601.35 


0.21 


0.20 


1.75 


0.20 


0.19 


11.00 


0.92 


0.74 


1 


620.00 


0.23 


0.26 


0.87 


0.16 


0.24 


5.90 


1.37 


1.01 


1 


641.11 


0.26 


0.25 


1.30 


0.22 


0.20 


8.39 


0.97 


0.96 


1 


669.51 


0.19 


0.20 


1.07 


0.21 


0.14 


9.43 


1.21 


0.70 


2 


359.35 


0.13 


0.15 


0.86 


0.15 


0.18 


6.15 


0.82 


1.11 


2 


390.57 


0.18 


0.24 


1.16 


0.31 


0.27 


4.84 


0.82 


0.54 


2 


422.54 


0.12 


0.13 


0.81 


0.10 


0.11 


7.95 


0.65 


0.69 


2 


454.72 


0.21 


0.22 


1.00 


0.12 


0.13 


10.85 


0.84 


1.07 


2 


486.64 


0.17 


0.16 


1.09 


0.22 


0.21 


7.16 


0.66 


1.42 


2 


518.58 


0.12 


0.12 


1.09 


0.14 


0.13 


9.85 


1.16 


0.68 


2 


551.08 


0.13 


0.14 


0.80 


0.13 


0.14 


10.34 


0.89 


0.96 


2 


584.11 


0.45 


0.59 


2.44 


0.21 


0.25 


8.50 


0.90 


0.65 


2 


616.71 


0.53 


0.38 


1.51 


0.25 


0.27 


8.24 


0.59 


0.62 


2 


648.92 


0.17 


0.18 


1.26 


0.17 


0.18 


7.06 


0.47 


0.62 
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